Near-unity Cooper pair splitting efficiency 
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The two electrons of a Cooper pair in a conventional superconductor form a singlet and therefore 
a maximally entangled state. Recently, it was demonstrated that the two particles can be extracted 
from the superconductor into two spatially separated contacts via two quantum dots (QDs) in a 
process called Cooper pair splitting (CPS). Competing transport processes, however, limit the effi- 
ciency of this process. Here we demonstrate efficiencies up to 90%, significantly larger than required 
to demonstrate interaction-dominated CPS, and on the right order to test Bell's inequality with 
electrons. We compare the CPS currents through both QDs, for which large apparent discrepan- 
cies are possible. The latter we explain intuitively and in a semi-classical master equation model. 
Large efficiencies are required to detect electron entanglement and for prospective electronics-based 
quantum information technologies. 
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Quantum entanglement between two particles is a fun- 
damental resource for quantum information technologies 
PP. Experiments with entangled photons are well devel- 
oped and already offer first applications [2J. However, 
entanglement between electrons, the fundamental parti- 
cles of electronics, is difficult to create in a controlled 
way. Electron-electron interactions, for example in the 
Fermi sea of a metal, tend to destroy particle correla- 
tions. In contrast, in nanostructured electronic devices 
interactions can be exploited, for example to generate 
entangled photons [3]. 

An elegant source of electronic entanglement are sin- 
glet two-particle ground states, for example in the nat- 
urally occurring BCS state of a conventional supercon- 
ductor. It was proposed to use such a superconductor to 
produce spatially separated entangled electrons in a pro- 
cess known as crossed Andreev reflection or Cooper pair 
splitting (CPS) [H[S]. Though metallic structures show 
electronic correlations due to superconductivity [BHlOj. 
their tunability is minimal. Recently, CPS was demon- 
strated on devices where a superconducting contact is 
coupled to two parallel quantum dots (QDs), each with 
a normal metal output lead, as shown schematically in 
Fig. la [XTl - TH] . In the latter experiments the CPS effi- 
ciency ranged from a few percent up to 50%. Such val- 
ues can in principle be reached without electron-electron 
interactions, e.g. in a chaotic cavity [15], or in a double- 
dot system with strong inter-dot coupling [T2] , where the 
electrons of a Cooper pair can exit the device through two 
ports at random. However, for applications and more so- 
phisticated experiments, for example the explicit demon- 
stration of entanglement, efficiencies close to unity are 
required. 

Here we present CPS experiments on a carbon nan- 
otube (CNT) device and demonstrate efficiencies up to 
90%, values only possible with electron-electron interac- 
tions. In addition, we find discrepancies when extracting 
the CPS part of the currents through the two quantum 



dots, which we relate to a competition between local 
processes and CPS in a semi-classical master equation 
model. The large CPS efficiencies and the increased un- 
derstanding of the relevant mechanisms are important 
steps on the way to an all-electronic source of entangled 
electron pairs in a solid-state device. 

An artificially colored scanning electron micrograph 
of a CPS device is shown in Fig. lb, together with a 
schematic of the measurement. A CVD-grown CNT (ar- 
row) is contacted in the center by an aluminum contact 
(S), which becomes superconducting below ~ 1.1 K and 
is evaporated on a 4nm palladium (Pd) contact layer. 
Two Pd contacts to the right and left of S serve as nor- 
mal metal contacts Nl and N2, both of which define a 
quantum dot (QD1 and QD2) on the two CNT segments 
adjacent to S. The QDs can be tuned electrically by a 
global backgate and the local side-gates SGI and SG2. 

From standard charge stability diagrams we extract 
charging energies of ~ 7meV for QD1 and ~ 4meV for 
QD2, an orbital energy spacing of ~ 1 meV, and an en- 
ergy gap due to the superconductor of ~ 120 fieV. With 
S in the normal state we find typical level broadenings 
of ~ 150-500 fieV. Relatively low peak conductances sug- 
gest rather asymmetric coupling of the QDs to the leads. 
The lever arms from a side-gate across the superconduc- 
tor to the other QD is roughly ten times smaller than that 




FIG. 1: (Color online) (a) Schematic of the device and rele- 
vant process probabilities, (b) Scanning electron micrograph 
of the CPS device and measurement setup. 
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FIG. 2: (a) Differential conductance of QD1, G\, and (b) of 
QD2, G2, as a function of the side-gate voltages Usgi and 
t^SG2- (c) and (d) Energy diagrams of local pair tunneling 
(LPT) and CPS. 



of a local side gate. In conductance measurements with 
the two QDs in series we do not observe an anti-crossing 
of the QD resonances, which suggests that the direct 
inter-dot tunnel coupling is negligible compared to pro- 
cesses via the superconductor. The experiments are per- 
formed in a dilution refrigerator at a base-temperature 
of - 20 mK. 

Figures 2a and 2b show the simultaneously recorded 
differential conductances G\ through QD1 and G 2 
through QD2, both as a function of the side-gate volt- 
ages t/sGl and t/sG2- The measurements were done at 
zero bias and zero magnetic field. When U^gi is varied, 
QD1 is tuned through several resonances, which result 
in conductance maxima in G\, labeled LI, L2 and L3 in 
Fig. 2a. The amplitudes of the resonances vary only little 
when tuning £7sG2, while the resonance position changes 
slightly due to capacitive cross talk from SG2 to QD1. 
Very weak, but similar conductance ridges labeled Rl, 
R2 and R3 can be observed in the conductance through 
QD2 in Fig. 2b. These are mainly tuned by SG2, which 
results in conductance ridges almost perpendicular to the 
ones in Fig. 2a due to QD1. 

Our main experimental findings are pronounced peaks 
when both QDs are in resonance. At these gate config- 
urations the conductance is increased by up to a factor 
of ~ 100 compared to the respective conductance ridge. 
This is most prominent in G2, but most of the peaks 
can also be observed in G\ on a larger background. No 
peaks at resonance crossings can be observed when the 
superconductivity is suppressed by a small external mag- 
netic field (see below). If only one QD is resonant, only 
local transport through this QD is allowed. A possible 



local process is local pair tunneling (LPT), illustrated in 
Fig. 2c: the first electron of a Cooper pair is emitted into 
the QD, which leaves S in a virtual excited state. When 
the first electron has left the dot, the second tunnels into 
the same QD. Other local processes like double charging 
of a dot are strongly suppressed by the large charging en- 
ergies. However, if both QDs are in resonance, the second 
electron can tunnel into QD2, as shown in Fig. 2d, which 
splits the initial Cooper pair. 

We now focus on the resonance crossing (L2,R2). Fig- 
ure 3a shows the Coulomb blockade resonance L2 in Gi 
as a function of J/sgi (bottom curve). In the same gate 
sweep, G2 is tuned through the resonance R2 due to ca- 
pacitive cross-talk, which results in a wide conductance 
maximum. However, an additional much sharper peak 
occurs at the voltage of the L2 resonance, with similar 
width and shape as the resonance in G\. When the su- 
perconductivity is suppressed by an external magnetic 
field of 250 mT, we find no additional peak in G2 at 
the resonance crossing, but a slight reduction consistent 
with a classical resistor network [TT] , see inset of Fig. 3a. 
The resonance positions do not change with field, but 
the overall conductance can vary strongly due to the su- 
perconductor's gap, which reduces local single electron 
transport and favors LPT. 

To assess CPS in the experiments we use the amplitude 
AG2 of the additional peak in G 2 at the position of the 
QD1 resonance, as illustrated in Fig. 3a. The subtracted 
background is determined by manually interpolating the 
bare QD2 resonance. Also indicated in Fig. 3a is the de- 
tuning SU between the two resonances. Figure 3b shows 
a series of SGl-sweeps at different values of U$G2 near the 
resonance crossing (L2,R2), with the curve from Fig. 3a 
highlighted in red. One finds that AG2 depends strongly 
on the detuning SU. In Fig. 3c we therefore plot AG2 vs. 
SU, where the value of Fig. 3a is marked by a red trian- 
gle. As another example, the conductance variation near 
the crossing (L3,R2) is also plotted in Fig. 3c. For all 
crossings we find that AG2 has a maximum at SU rj 0, 
i.e. where both QDs arc in resonance, in agreement with 
theoretical predictions [4]. For SU =^ 0, AG2 falls off 
rapidly and tends to zero on an energy scale consistent 
with the width of the respective resonances. 

On the resonance crossings investigated here the max- 
imum change in G2 is 0.012 e 2 /h. This number has 
to be compared to the total conductance, including the 
local processes, so that we define the visibility of CPS 
in the second branch of the Cooper pair splitter as 
772 = AG2/G2 (similar for G x ). The CPS visibilities for 
both branches on resonance crossing (L3,R2) are plotted 
in Fig. 3d. r\i is essentially constant over a large range 
of SU and reaches values up to 98%, i.e. the current in 
one branch can be dominated by CPS. rjx, however, has 
a maximum of only 73% at SU « and drops to zero for 
a large detuning. 

As a measure for the CPS efficiency we compare the 
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FIG. 3: (Color online) (a) Gi and G2 as a function of (7sgi 
for UsG2 ~ 2.07 V. Inset: the same measurement at a mag- 
netic field of 250 mT to suppress the superconductivity, (b) 
G2(Usgi) for a series of side-gate voltages Usg2- (c) AG2 as 
a function of the detuning SU between the resonances in QD1 
and QD2 for the indicated resonance-crossings, (d) Plots of 
the visibilities rji, the CPS efficiency s and rj = 771772. 



CPS currents to the total currents in both branches of 
the device. Assuming that CPS leads to a conductance 
Gcps m each branch, independent of other processes, we 
define the CPS efficiency as 



2G 



CPS 



G\ + G2 



(1) 



By setting Gcps = AG2 we find efficiencies up to s rj 
90%, much larger than required to demonstrate interac- 
tion dominated CPS. The efficiency as a function of SU 
is plotted in Fig. 3d for the crossing (L3,R2). However, 
depending on the intended purpose of the entangler, s is 
not necessarily the relevant parameter. For example, in 
tests of Bell's inequality proposed for electrons [TBI [T7] . 
the measured quantities are current cross correlations be- 
tween the normal metal terminals, which suggests to use 
the following figure of merit: 



V = Vi ■ % 



AG 1 AG 2 
Gi G2 



(2) 



A violation of Bell's inequality requires 77 > l/v2 « 71%. 
In Fig. 3d, 77 is plotted as a function of SU for the cross- 
ing (L3,R2). We find values up to 77 = 68%, mostly lim- 
ited by the large rates of local processes through QD1. 
Nontheless, the large visibility in G2 demonstrates the 



feasibility of testing Bell's inequality with electrons, if an 
ideal detection scheme was available. 

Intuitively one might expect AGi = AG2. This is 
found within experimental errors for 4 of the 9 resonance 
crossings. As an example, AGi and AG2 of the crossing 
(L3,R2) investigated above are plotted as a function of 
SU in Fig. 4a. For the other crossings, the two conduc- 
tance variations deviate significantly from each other. 4 
of the 9 crossings exhibit curves comparable to (L3,R1) 
plotted in Fig. 4b. Here, AG2 is larger than AGi by 
about a factor of 2, but with a similar curve shape. One 
of the 9 crossings, (L2,R1) shown in Fig. 4c, is very pe- 
culiar: the variation in AGi is almost negligible, while 
AG2 exhibits a pronounced peak. In addition, one finds 
that AG2 > Gi, i.e. the conductance variation in one 
branch is larger than the total conductance in the other. 

To explain our experiments we discuss a strongly sim- 
plified semi-classical master equation model. More so- 
phisticated models can be found in [19l |20]. For each 
QD we consider a single level with constant broadening 
and a large charging energy. The system can be in one 
of the following four states: both QDs empty, either QD 
filled with one electron, or both dots occupied. The tun- 
neling processes illustrated in Fig. la lead to transitions 
between these states. We assume that effectively elec- 
trons are transfered only in one direction, from S to the 
QDs and from the QDs to the respective normal metal 
contact. In addition we consider a tunnel coupling be- 
tween the dots. The pi in Fig. la are the probabilities 
that the corresponding process changes the occupation of 
a system state. It is not trivial to extract absolute values 
for the pi from the experiments, especially for the com- 
plex transport processes involving S. The resonances are 
incorporated as a normalized effective density of states. 
We use a diagrammatic method based on maximal trees 
to obtain the steady-state occupation probabilities from 
the corresponding master equation |18j . From the popu- 
lations of the QDs we then calculate the transport rates. 

Our model reproduces qualitatively the observed con- 
ductance variations and shows that a finite QD popu- 
lation can lead to a competition between the various 
transport mechanisms. In Figs. 4d-f simulated conduc- 
tance variations are plotted for different QD1 parameters, 
while QD2 is kept at p§2 = 0.01 and j»N2 = 0.1, i.e. in 
the regime of Ref. [J], where the coupling to S is much 
weaker than to the normal contacts. We set pcps = 0.03 
to obtain relative amplitudes comparable to the experi- 
ments, and P12 = 0.001 so that the inter-dot coupling is 
the smallest parameter in the problem. If the two cur- 
rent branches are similar, i.e. psi ~ ps2 < Pni ~ P~N2, 
one finds AGi = AG2, as shown in Fig. 4d [21] . similar 
to the experiments presented in Fig. 4a. 

For asymmetric branches, the conductance variations 
are not identical anymore. Figure 4e shows plots for 
Psi = Pni = 0.1, for which AG2 ~ 2AGi, as in the 
experiment shown in Fig. 4b. The model also allows us 
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FIG. 4: (Color online) (a)-(c) AGi and AG2 as a function of 
the detuning 8U for the indicated resonance crossings, (d)- 
(f) Similar plots obtained from the master equation model 
of CPS, including Gcps- The parameters varied between the 
simulations are given in the figures. The inset in (e) shows 
the conductances in the branch of QD1 due to local processes 
(Gsi), CPS and at Nl. 



to calculate the rate at which Cooper pairs are extracted 
from S by CPS. The corresponding conductance is also 
plotted in Figs. 4d-f. We find that AG, < Gcps as long 
as the inter-dot coupling p\2 is negligible, i.e. the experi- 
mentally extracted CPS conductance underestimates the 
actual value, max(AGi, AG2) < Gcps- The explanation 
is that due to CPS on a resonance crossing the average 
QD populations are increased beyond the off-resonance 
equilibrium due to the local processes, which leads to a 
reduction of the current into the N contacts. This is illus- 
trated in the inset of Fig. 4e, where the calculated local 
conductance from SI, Gsi, has a minimum where Gcps 
is maximal. Intuitively, the QDs are not emptied fast 
enough, which inhibits all processes on the dot. 

The situation is more complex if the tunnel coupling 
between the dots becomes relevant. For example, if 
Pni = P12 — 0.001 and psi = 0.01, as used in Fig. 4f, the 
electrons can leave QD1 to Nl and to QD2 with the same 
probability. Since pni is small, this quenches Gi, but 
G2 is increased due to the additional current from QD1. 
Here, the AG^ do not give an upper or lower bound for 
the CPS rate and AG2 can become larger than Gi, as in 
the experimental curves in Fig. 4c. We note that the dis- 
cussed situations are not in the regime of completely filled 
QDs. Our model suggests that in this unitary limit the 
conductances can be reduced considerably in the center 
of a resonance crossing. In our data we find no evidence 
for this prediction. However, since the coupling between 
S and an InAs nanowires can be strong, the model might 



account for the as yet unexplained anomalous behavior 
of the on-resonance signals in [TTJ . 

In summary, we present Cooper pair splitting exper- 
iments with efficiencies up to 90%, demonstrating the 
importance of electron-electron interactions in such sys- 
tems. For the figure of merit relevant in tests of Bell's 
inequality for electons we find values close to the required 
limits. In addition, we asses CPS on both QDs and find 
rather large apparent discrepancies between the two con- 
ductance variations, which we explain in a semi-classical 
master equation model. The latter suggests that for neg- 
ligible inter-dot couplings the experimentally extracted 
CPS rates are a lower bound to the real CPS rates. Our 
experiments and calculations show that there is a large 
variety of different transport phenomena in a Cooper pair 
splitting device that need further investigation. Of cap- 
ital importance is the observation that if both dots had 
the properties of QD2, tests of Bell's inequality even with 
non-ideal detectors could be performed to detect electron 
entanglement, an imprtant step on the way to a source 
of entangled electron pairs on demand. 
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